Insight into the adaptive role of arachnid genome-wide duplication through chromosome-level genome assembly of the Western black widow spider.

Although spiders are one of the most diverse groups of arthropods, the genetic architecture of their evolutionary adaptations is largely unknown. Specifically, ancient genome-wide duplication occurring during arachnid evolution ~450 mya resulted in a vast assembly of gene families, yet the extent to which selection has shaped this variation is understudied. To aid in comparative genome sequence analyses, we provide a chromosome-level genome of the Western black widow spider (Latrodectus hesperus)-a focus due to its silk properties, venom applications, and as a model for urban adaptation. We used long-read and Hi-C sequencing data, combined with transcriptomes, to assemble 14 chromosomes in a 1.46 Gb genome, with 38,393 genes annotated, and a BUSCO score of 95.3%. Our analyses identified high repetitive gene content and heterozygosity, consistent with other spider genomes, which has led to challenges in genome characterization. Our comparative evolutionary analyses of eight genomes available for species within the Araneoidea group (orb weavers and their descendants) identified 1,827 single-copy orthologs. Of these, 155 exhibit significant positive selection primarily associated with developmental genes, and with traits linked to sensory perception. These results support the hypothesis that several traits unique to spiders emerged from the adaptive evolution of ohnologs-or retained ancestrally duplicated genes-from ancient genome-wide duplication. These comparative spider genome analyses can serve as a model to understand how positive selection continually shapes ancestral duplications in generating novel traits today within and between diverse taxonomic groups.

Correlation between protein secondary structure and mechanical performance for the ultra-tough dragline silk of Darwin's bark spider.

The spider major ampullate (MA) silk exhibits high tensile strength and extensibility and is typically a blend of MaSp1 and MaSp2 proteins with the latter comprising glycine-proline-glycine-glycine-X repeating motifs that promote extensibility and supercontraction. The MA silk from Darwin's bark spider (Caerostris darwini) is estimated to be two to three times tougher than the MA silk from other spider species. Previous research suggests that a unique MaSp4 protein incorporates proline into a novel glycine-proline-glycine-proline motif and may explain C. darwini MA silk's extraordinary toughness. However, no direct correlation has been made between the silk's molecular structure and its mechanical properties for C. darwini. Here, we correlate the relative protein secondary structure composition of MA silk from C. darwini and four other spider species with mechanical properties before and after supercontraction to understand the effect of the additional MaSp4 protein. Our results demonstrate that C. darwini MA silk possesses a unique protein composition with a lower ratio of helices (31%) and ß-sheets (20%) than other species. Before supercontraction, toughness, modulus and tensile strength correlate with percentages of ß-sheets, unordered or random coiled regions and ß-turns. However, after supercontraction, only modulus and strain at break correlate with percentages of ß-sheets and ß-turns. Our study highlights that additional information including crystal size and crystal and chain orientation is necessary to build a complete structure-property correlation model.

Expanding Canonical Spider Silk Properties through a DNA Combinatorial Approach.

The properties of native spider silk vary within and across species due to the presence of different genes containing conserved repetitive core domains encoding a variety of silk proteins. Previous studies seeking to understand the function and material properties of these domains focused primarily on the analysis of dragline silk proteins, MaSp1 and MaSp2. Our work seeks to broaden the mechanical properties of silk-based biomaterials by establishing two libraries containing genes from the repetitive core region of the native Latrodectus hesperus silk genome (Library A: genes masp1, masp2, tusp1, acsp1; Library B: genes acsp1, pysp1, misp1, flag). The expressed and purified proteins were analyzed through Fourier Transform Infrared Spectrometry (FTIR). Some of these new proteins revealed a higher portion of ß-sheet content in recombinant proteins produced from gene constructs containing a combination of masp1/masp2 and acsp1/tusp1 genes than recombinant proteins which consisted solely of dragline silk genes (Library A). A higher portion of ß-turn and random coil content was identified in recombinant proteins from pysp1 and flag genes (Library B). Mechanical characterization of selected proteins purified from Library A and Library B formed into films was assessed by Atomic Force Microscopy (AFM) and suggested Library A recombinant proteins had higher elastic moduli when compared to Library B recombinant proteins. Both libraries had higher elastic moduli when compared to native spider silk proteins. The preliminary approach demonstrated here suggests that repetitive core regions of the aforementioned genes can be used as building blocks for new silk-based biomaterials with varying mechanical properties.

Spider phylosymbiosis: divergence of widow spider species and their tissues' microbiomes.

BACKGROUND: Microbiomes can have profound impacts on host biology and evolution, but to date, remain vastly understudied in spiders despite their unique and diverse predatory adaptations. This study evaluates closely related species of spiders and their host-microbe relationships in the context of phylosymbiosis, an eco-evolutionary pattern where the microbial community profile parallels the phylogeny of closely related host species. Using 16S rRNA gene amplicon sequencing, we characterized the microbiomes of five species with known phylogenetic relationships from the family Theridiidae, including multiple closely related widow spiders (L. hesperus, L. mactans, L. geometricus, S. grossa, and P. tepidariorum). RESULTS: We compared whole animal and tissue-specific microbiomes (cephalothorax, fat bodies, venom glands, silk glands, and ovary) in the five species to better understand the relationship between spiders and their microbial symbionts. This showed a strong congruence of the microbiome beta-diversity of the whole spiders, cephalothorax, venom glands, and silk glands when compared to their host phylogeny. Our results support phylosymbiosis in these species and across their specialized tissues. The ovary tissue microbial dendrograms also parallel the widow phylogeny, suggesting vertical transfer of species-specific bacterial symbionts. By cross-validating with RNA sequencing data obtained from the venom glands, silk glands and ovaries of L. hesperus, L. geometricus, S. grossa, and P. tepidariorum we confirmed that several microbial symbionts of interest are viably active in the host. CONCLUSION: Together these results provide evidence that supports the importance of host-microbe interactions and the significant role microbial communities may play in the evolution and adaptation of their hosts.

Structural and functional characterization of sticholysin III: A newly discovered actinoporin within the venom of the sea anemone Stichodactyla helianthus.

Actinoporins are a family of pore-forming toxins produced by sea anemones as part of their venomous cocktail. These proteins remain soluble and stably folded in aqueous solution, but when interacting with sphingomyelin-containing lipid membranes, they become integral oligomeric membrane structures that form a pore permeable to cations, which leads to cell death by osmotic shock. Actinoporins appear as multigenic families within the genome of sea anemones: several genes encoding very similar actinoporins are detected within the same species. The Caribbean Sea anemone Stichodactyla helianthus produces three actinoporins (sticholysins I, II and III; StnI, StnII and StnIII) that differ in their toxic potency. For example, StnII is about four-fold more effective than StnI against sheep erythrocytes in causing hemolysis, and both show synergy. However, StnIII, recently discovered in the S. helianthus transcriptome, has not been characterized so far. Here we describe StnIII's spectroscopic and functional properties and show its potential to interact with the other Stns. StnIII seems to maintain the well-preserved fold of all actinoporins, characterized by a high content of ß-sheet, but it is significantly less thermostable. Its functional characterization shows that the critical concentration needed to form active pores is higher than for either StnI or StnII, suggesting differences in behavior when oligomerizing on membrane surfaces. Our results show that StnIII is an interesting and unexpected piece in the puzzle of how this Caribbean Sea anemone species modulates its venomous activity.

Ovarian Transcriptomic Analyses in the Urban Human Health Pest, the Western Black Widow Spider.

Due to their abundance and ability to invade diverse environments, many arthropods have become pests of economic and health concern, especially in urban areas. Transcriptomic analyses of arthropod ovaries have provided insight into life history variation and fecundity, yet there are few studies in spiders despite their diversity within arthropods. Here, we generated a de novo ovarian transcriptome from 10 individuals of the western black widow spider (Latrodectus hesperus), a human health pest of high abundance in urban areas, to conduct comparative ovarian transcriptomic analyses. Biological processes enriched for metabolism-specifically purine, and thiamine metabolic pathways linked to oocyte development-were significantly abundant in L. hesperus. Functional and pathway annotations revealed overlap among diverse arachnid ovarian transcriptomes for highly-conserved genes and those linked to fecundity, such as oocyte maturation in vitellogenin and vitelline membrane outer layer proteins, hormones, and hormone receptors required for ovary development, and regulation of fertility-related genes. Comparative studies across arachnids are greatly needed to understand the evolutionary similarities of the spider ovary, and here, the identification of ovarian proteins in L. hesperus provides potential for understanding how increased fecundity is linked to the success of this urban pest.

Alternative Transcription at Venom Genes and Its Role as a Complementary Mechanism for the Generation of Venom Complexity in the Common House Spider.

The complex composition of venom, a proteinaceous secretion used by diverse animal groups for predation or defense, is typically viewed as being driven by gene duplication in conjunction with positive selection, leading to large families of diversified toxins with selective venom gland expression. Yet, the production of alternative transcripts at venom genes is often overlooked as another potentially important process that could contribute proteins to venom, and requires comprehensive datasets integrating genome and transcriptome sequences together with proteomic characterization of venom to be fully documented. In the common house spider, Parasteatoda tepidariorum, we used RNA sequencing of four tissue types in conjunction with the sequenced genome to provide a comprehensive transcriptome annotation. We also used mass spectrometry to identify a minimum of 99 distinct proteins in P tepidariorum venom, including at least 33 latrotoxins, pore-forming neurotoxins shared with the confamilial black widow. We found that venom proteins are much more likely to come from multiple transcript genes, whose transcripts produced distinct protein sequences. The presence of multiple distinct proteins in venom from transcripts at individual genes was confirmed for eight loci by mass spectrometry, and is possible at 21 others. Alternative transcripts from the same gene, whether encoding or not encoding a protein found in venom, showed a range of expression patterns, but were not necessarily restricted to the venom gland. However, approximately half of venom protein encoding transcripts were found among the 1,318 transcripts with strongly venom gland biased expression. Our findings revealed an important role for alternative transcription in generating venom protein complexity and expanded the traditional model of venom evolution.

The transcriptome of Darwin's bark spider silk glands predicts proteins contributing to dragline silk toughness.

Darwin's bark spider (Caerostris darwini) produces giant orb webs from dragline silk that can be twice as tough as other silks, making it the toughest biological material. This extreme toughness comes from increased extensibility relative to other draglines. We show C. darwini dragline-producing major ampullate (MA) glands highly express a novel silk gene transcript (MaSp4) encoding a protein that diverges markedly from closely related proteins and contains abundant proline, known to confer silk extensibility, in a unique GPGPQ amino acid motif. This suggests C. darwini evolved distinct proteins that may have increased its dragline's toughness, enabling giant webs. Caerostris darwini's MA spinning ducts also appear unusually long, potentially facilitating alignment of silk proteins into extremely tough fibers. Thus, a suite of novel traits from the level of genes to spinning physiology to silk biomechanics are associated with the unique ecology of Darwin's bark spider, presenting innovative designs for engineering biomaterials.

Stichodactyla helianthus' de novo transcriptome assembly: Discovery of a new actinoporin isoform.

Transcriptomic profiling of venom producing tissues from different animals is an effective approach for discovering new toxins useful in biotechnological and pharmaceutical applications, as well in evolutionary comparative studies of venomous animals. Stichodactyla helianthus is a Caribbean sea anemone which produces actinoporins as part of its toxic venom. This family of pore forming toxins is multigenic and at least two different isoforms, encoded by separate genes, are produced by S. helianthus. These isoforms, sticholysins I and II, share 93% amino acid identity but differ in their pore forming activity and act synergistically. This observation suggests that other actinoporin isoforms, if present in the venomous mixture, could offer an advantageous strategy to modulate whole venom activity. Using high-throughput sequencing we generated a de novo transcriptome of S. helianthus and determined the relative expression of assembled transcripts using RNA-Seq to better characterize components of this species' venom, focusing on actinoporin diversity. Applying this approach, we have discovered at least one new actinoporin variant from S. helianthus in addition to several other putative venom components.

Recent progress and prospects for advancing arachnid genomics.

Arachnids exhibit tremendous species richness and adaptations of biomedical, industrial, and agricultural importance. Yet genomic resources for arachnids are limited, with the first few spider and scorpion genomes becoming accessible in the last four years. We review key insights from these genome projects, and recommend additional genomes for sequencing, emphasizing taxa of greatest value to the scientific community. We suggest greater sampling of spiders whose genomes are understudied but hold important protein recipes for silk and venom production. We further recommend arachnid genomes to address significant evolutionary topics, including the phenotypic impact of genome duplications. A barrier to high-quality arachnid genomes are assemblies based solely on short-read data, which may be overcome by long-range sequencing and other emerging methods.

Evolutionary shifts in gene expression decoupled from gene duplication across functionally distinct spider silk glands.

Spider silk synthesis is an emerging model for the evolution of tissue-specific gene expression and the role of gene duplication in functional novelty, but its potential has not been fully realized. Accordingly, we quantified transcript (mRNA) abundance in seven silk gland types and three non-silk gland tissues for three cobweb-weaving spider species. Evolutionary analyses based on expression levels of thousands of homologous transcripts and phylogenetic reconstruction of 605 gene families demonstrated conservation of expression for each gland type among species. Despite serial homology of all silk glands, the expression profiles of the glue-forming aggregate glands were divergent from fiber-forming glands. Also surprising was our finding that shifts in gene expression among silk gland types were not necessarily coupled with gene duplication, even though silk-specific genes belong to multi-paralog gene families. Our results challenge widely accepted models of tissue specialization and significantly advance efforts to replicate silk-based high-performance biomaterials.

Duplication and concerted evolution of MiSp-encoding genes underlie the material properties of minor ampullate silks of cobweb weaving spiders.

BACKGROUND: Orb-web weaving spiders and their relatives use multiple types of task-specific silks. The majority of spider silk studies have focused on the ultra-tough dragline silk synthesized in major ampullate glands, but other silk types have impressive material properties. For instance, minor ampullate silks of orb-web weaving spiders are as tough as draglines, due to their higher extensibility despite lower strength. Differences in material properties between silk types result from differences in their component proteins, particularly members of the spidroin (spider fibroin) gene family. However, the extent to which variation in material properties within a single silk type can be explained by variation in spidroin sequences is unknown. Here, we compare the minor ampullate spidroins (MiSp) of orb-weavers and cobweb weavers. Orb-web weavers use minor ampullate silk to form the auxiliary spiral of the orb-web while cobweb weavers use it to wrap prey, suggesting that selection pressures on minor ampullate spidroins (MiSp) may differ between the two groups. RESULTS: We report complete or nearly complete MiSp sequences from five cobweb weaving spider species and measure material properties of minor ampullate silks in a subset of these species. We also compare MiSp sequences and silk properties of our cobweb weavers to published data for orb-web weavers. We demonstrate that all our cobweb weavers possess multiple MiSp loci and that one locus is more highly expressed in at least two species. We also find that the proportion of ß-spiral-forming amino acid motifs in MiSp positively correlates with minor ampullate silk extensibility across orb-web and cobweb weavers. CONCLUSIONS: MiSp sequences vary dramatically within and among spider species, and have likely been subject to multiple rounds of gene duplication and concerted evolution, which have contributed to the diverse material properties of minor ampullate silks. Our sequences also provide templates for recombinant silk proteins with tailored properties.

House spider genome uncovers evolutionary shifts in the diversity and expression of black widow venom proteins associated with extreme toxicity.

BACKGROUND: Black widow spiders are infamous for their neurotoxic venom, which can cause extreme and long-lasting pain. This unusual venom is dominated by latrotoxins and latrodectins, two protein families virtually unknown outside of the black widow genus Latrodectus, that are difficult to study given the paucity of spider genomes. Using tissue-, sex- and stage-specific expression data, we analyzed the recently sequenced genome of the house spider (Parasteatoda tepidariorum), a close relative of black widows, to investigate latrotoxin and latrodectin diversity, expression and evolution. RESULTS: We discovered at least 47 latrotoxin genes in the house spider genome, many of which are tandem-arrayed. Latrotoxins vary extensively in predicted structural domains and expression, implying their significant functional diversification. Phylogenetic analyses show latrotoxins have substantially duplicated after the Latrodectus/Parasteatoda split and that they are also related to proteins found in endosymbiotic bacteria. Latrodectin genes are less numerous than latrotoxins, but analyses show their recruitment for venom function from neuropeptide hormone genes following duplication, inversion and domain truncation. While latrodectins and other peptides are highly expressed in house spider and black widow venom glands, latrotoxins account for a far smaller percentage of house spider venom gland expression. CONCLUSIONS: The house spider genome sequence provides novel insights into the evolution of venom toxins once considered unique to black widows. Our results greatly expand the size of the latrotoxin gene family, reinforce its narrow phylogenetic distribution, and provide additional evidence for the lateral transfer of latrotoxins between spiders and bacterial endosymbionts. Moreover, we strengthen the evidence for the evolution of latrodectin venom genes from the ecdysozoan Ion Transport Peptide (ITP)/Crustacean Hyperglycemic Hormone (CHH) neuropeptide superfamily. The lower expression of latrotoxins in house spiders relative to black widows, along with the absence of a vertebrate-targeting α-latrotoxin gene in the house spider genome, may account for the extreme potency of black widow venom.

Effects of Gene Duplication, Positive Selection, and Shifts in Gene Expression on the Evolution of the Venom Gland Transcriptome in Widow Spiders.

Gene duplication and positive selection can be important determinants of the evolution of venom, a protein-rich secretion used in prey capture and defense. In a typical model of venom evolution, gene duplicates switch to venom gland expression and change function under the action of positive selection, which together with further duplication produces large gene families encoding diverse toxins. Although these processes have been demonstrated for individual toxin families, high-throughput multitissue sequencing of closely related venomous species can provide insights into evolutionary dynamics at the scale of the entire venom gland transcriptome. By assembling and analyzing multitissue transcriptomes from the Western black widow spider and two closely related species with distinct venom toxicity phenotypes, we do not find that gene duplication and duplicate retention is greater in gene families with venom gland biased expression in comparison with broadly expressed families. Positive selection has acted on some venom toxin families, but does not appear to be in excess for families with venom gland biased expression. Moreover, we find 309 distinct gene families that have single transcripts with venom gland biased expression, suggesting that the switching of genes to venom gland expression in numerous unrelated gene families has been a dominant mode of evolution. We also find ample variation in protein sequences of venom gland-specific transcripts, lineage-specific family sizes, and ortholog expression among species. This variation might contribute to the variable venom toxicity of these species.

Spider Transcriptomes Identify Ancient Large-Scale Gene Duplication Event Potentially Important in Silk Gland Evolution.

The evolution of specialized tissues with novel functions, such as the silk synthesizing glands in spiders, is likely an influential driver of adaptive success. Large-scale gene duplication events and subsequent paralog divergence are thought to be required for generating evolutionary novelty. Such an event has been proposed for spiders, but not tested. We de novo assembled transcriptomes from three cobweb weaving spider species. Based on phylogenetic analyses of gene families with representatives from each of the three species, we found numerous duplication events indicative of a whole genome or segmental duplication. We estimated the age of the gene duplications relative to several speciation events within spiders and arachnids and found that the duplications likely occurred after the divergence of scorpions (order Scorpionida) and spiders (order Araneae), but before the divergence of the spider suborders Mygalomorphae and Araneomorphae, near the evolutionary origin of spider silk glands. Transcripts that are expressed exclusively or primarily within black widow silk glands are more likely to have a paralog descended from the ancient duplication event and have elevated amino acid replacement rates compared with other transcripts. Thus, an ancient large-scale gene duplication event within the spider lineage was likely an important source of molecular novelty during the evolution of silk gland-specific expression. This duplication event may have provided genetic material for subsequent silk gland diversification in the true spiders (Araneomorphae).

Extraction of venom and venom gland microdissections from spiders for proteomic and transcriptomic analyses.

Venoms are chemically complex secretions typically comprising numerous proteins and peptides with varied physiological activities. Functional characterization of venom proteins has important biomedical applications, including the identification of drug leads or probes for cellular receptors. Spiders are the most species rich clade of venomous organisms, but the venoms of only a few species are well-understood, in part due to the difficulty associated with collecting minute quantities of venom from small animals. This paper presents a protocol for the collection of venom from spiders using electrical stimulation, demonstrating the procedure on the Western black widow (Latrodectus hesperus). The collected venom is useful for varied downstream analyses including direct protein identification via mass spectrometry, functional assays, and stimulation of venom gene expression for transcriptomic studies. This technique has the advantage over protocols that isolate venom from whole gland homogenates, which do not separate genuine venom components from cellular proteins that are not secreted as part of the venom. Representative results demonstrate the detection of known venom peptides from the collected sample using mass spectrometry. The venom collection procedure is followed by a protocol for dissecting spider venom glands, with results demonstrating that this leads to the characterization of venom-expressed proteins and peptides at the sequence level.

Diverse formulas for spider dragline fibers demonstrated by molecular and mechanical characterization of spitting spider silk.

Spider silks have outstanding mechanical properties. Most research has focused on dragline silk proteins (major ampullate spidroins, MaSps) from orb-weaving spiders. Using silk gland expression libraries from the haplogyne spider Scytodes thoracica, we discovered two novel spidroins (S. thoracica fibroin 1 and 2). The amino acid composition of S. thoracica silk glands and dragline fibers suggest that fibroin 1 is the major component of S. thoracica dragline silk. Fibroin 1 is dominated by glycine-alanine motifs, and lacks sequence motifs associated with orb-weaver MaSps. We hypothesize fibroin 2 is a piriform or aciniform silk protein, based on amino acid composition, spigot morphology, and phylogenetic analyses. S. thoracica's dragline silk is less tough than previously reported, but is still comparable to other dragline silks. Our analyses suggest that dragline silk proteins evolved multiple times. This demonstrates that spider dragline silk is more diverse than previously understood, providing alternative high performance silk designs.

Gene structure, regulatory control, and evolution of black widow venom latrotoxins.

Black widow venom contains α-latrotoxin, infamous for causing intense pain. Combining 33 kb of Latrodectus hesperus genomic DNA with RNA-Seq, we characterized the α-latrotoxin gene and discovered a paralog, 4.5 kb downstream. Both paralogs exhibit venom gland specific transcription, and may be regulated post-transcriptionally via musashi-like proteins. A 4 kb intron interrupts the α-latrotoxin coding sequence, while a 10 kb intron in the 3' UTR of the paralog may cause non-sense-mediated decay. Phylogenetic analysis confirms these divergent latrotoxins diversified through recent tandem gene duplications. Thus, latrotoxin genes have more complex structures, regulatory controls, and sequence diversity than previously proposed.

Dramatic expansion of the black widow toxin arsenal uncovered by multi-tissue transcriptomics and venom proteomics.

BACKGROUND: Animal venoms attract enormous interest given their potential for pharmacological discovery and understanding the evolution of natural chemistries. Next-generation transcriptomics and proteomics provide unparalleled, but underexploited, capabilities for venom characterization. We combined multi-tissue RNA-Seq with mass spectrometry and bioinformatic analyses to determine venom gland specific transcripts and venom proteins from the Western black widow spider (Latrodectus hesperus) and investigated their evolution. RESULTS: We estimated expression of 97,217 L. hesperus transcripts in venom glands relative to silk and cephalothorax tissues. We identified 695 venom gland specific transcripts (VSTs), many of which BLAST and GO term analyses indicate may function as toxins or their delivery agents. ~38% of VSTs had BLAST hits, including latrotoxins, inhibitor cystine knot toxins, CRISPs, hyaluronidases, chitinase, and proteases, and 59% of VSTs had predicted protein domains. Latrotoxins are venom toxins that cause massive neurotransmitter release from vertebrate or invertebrate neurons. We discovered ≥ 20 divergent latrotoxin paralogs expressed in L. hesperus venom glands, significantly increasing this biomedically important family. Mass spectrometry of L. hesperus venom identified 49 proteins from VSTs, 24 of which BLAST to toxins. Phylogenetic analyses showed venom gland specific gene family expansions and shifts in tissue expression. CONCLUSIONS: Quantitative expression analyses comparing multiple tissues are necessary to identify venom gland specific transcripts. We present a black widow venom specific exome that uncovers a trove of diverse toxins and associated proteins, suggesting a dynamic evolutionary history. This justifies a reevaluation of the functional activities of black widow venom in light of its emerging complexity.

Multi-tissue transcriptomics of the black widow spider reveals expansions, co-options, and functional processes of the silk gland gene toolkit.

BACKGROUND: Spiders (Order Araneae) are essential predators in every terrestrial ecosystem largely because they have evolved potent arsenals of silk and venom. Spider silks are high performance materials made almost entirely of proteins, and thus represent an ideal system for investigating genome level evolution of novel protein functions. However, genomic level resources remain limited for spiders. RESULTS: We de novo assembled a transcriptome for the Western black widow (Latrodectus hesperus) from deeply sequenced cDNAs of three tissue types. Our multi-tissue assembly contained ~100,000 unique transcripts, of which > 27,000 were annotated by homology. Comparing transcript abundance among the different tissues, we identified 647 silk gland-specific transcripts, including the few known silk fiber components (e.g. six spider fibroins, spidroins). Silk gland specific transcripts are enriched compared to the entire transcriptome in several functions, including protein degradation, inhibition of protein degradation, and oxidation-reduction. Phylogenetic analyses of 37 gene families containing silk gland specific transcripts demonstrated novel gene expansions within silk glands, and multiple co-options of silk specific expression from paralogs expressed in other tissues. CONCLUSIONS: We propose a transcriptional program for the silk glands that involves regulating gland specific synthesis of silk fiber and glue components followed by protecting and processing these components into functional fibers and glues. Our black widow silk gland gene repertoire provides extensive expansion of resources for biomimetic applications of silk in industry and medicine. Furthermore, our multi-tissue transcriptome facilitates evolutionary analysis of arachnid genomes and adaptive protein systems.